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Abstract. 

o: 

CN . We solve the equivalence problem for the orthogonally separable webs 

Ph! on the three-sphere under the action of the isometry group. This con- 

tinues a classical project initiated by Olevsky in which he solved the 
corresponding canonical forms problem. The solution to the equivalence 
problem together with the results by Olevsky forms a complete solution 
to the problem of orthogonal separation of variables to the Hamilton- 
& '. Jacobi equation defined on the three-sphere via orthogonal separation of 

variables. It is based on invariant properties of the characteristic Killing 
^ ; two-tensors in addition to properties of the corresponding algebraic cur- 

vature tensor and the associated Ricci tensor. The result is illustrated by 
a non-trivial application to a natural Hamiltonian defined on the three- 
sphere. 
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This paper presents a complete solution to the equivalence problem for orthogonally 
separable webs (OSWs) generated by the characteristic Killing tensors (CKTs) de- 
fined on the three-dimensional sphere § 3 . In addition, we lay the ground work for 
solving the analogous problem for three-dimensional hyperbolic space H 3 as well as 
the spaces of constant non-zero curvature of higher dimensions. This work continues 
^ ■ and generalizes our studies of the OSWs defined in flat (pseudo-)Riemannian spaces 

(see, for instance, [TJ [HJ [TB] and the relevant references therein). 

Reportedly, the study of OSWs generated by Killing two-tensors was conceived 
in various classical articles published throughout the 19th century. A list of such 
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works should include, among many others, Liouville [21], Neumann [26], Bertrand 
[3], Morera [25] . and Stackel [31]. The research in the area continued throughout the 
last century to make this study more systematic and develop new results. Of note 
was the celebrated 1934 paper by Eisenhart [8]. In this paper, the author completely 
solved the canonical forms problem for the OSWs generated by CKTs defined in E 3 by 
showing that there were exactly eleven inequivalent (in some sense) such OSWs and 
represented each of them with the corresponding canonical form. Furthermore, he 
gave a general criterion linking orthogonal separability with geometric and algebraic 
properties of Killing two-tensors. In addition to E 3 , he also partially solved the 
canonical forms problem for S 3 . This problem, as well as the case of hyperbolic space 
H 3 , was later completed by Olevsky in 1950. In his paper [28], Olevsky determined 
the number of orbits corresponding to the orthogonal separable coordinates on § 3 
and H 3 , listing the metrics in these coordinates and the transformations to separable 
coordinates in each case. Kalnins et al [20] completed Olevsky's work on S 3 by also 
listing the CKTs (using the language of differential operators) in their respective 
canonical forms representing the OSWs. In recent years, the research in the area has 
been successfully continued yielding many important results in the area in connection 
with the study of integrable and superintegrable classical and quantum Hamiltonian 
systems (see, for example, [U [HI [161 EH USB [35] and the relevant references therein). 

It must be noted, however, that many of the aforementioned investigations con- 
cerned the canonical forms problems for the CKTs in question, while the applications 
arising in the field of classical and quantum Hamiltonian systems often involve the 
properties of the CKTs in their general form as far as the action of the corresponding 
isometry groups are concerned (see [H], for example). From this perspective, the 
main goal of this paper is to extend the results by Olevsky and others and solve the 
equivalence problem for the CKTs defined on § 3 . 

The paper is organized as follows. Section 2 contains a brief review of the required 
mathematical tools defined in the framework of the invariant theory of Killing tensors 
(ITKTs) that will be used in solving the main problem. In Section 3 we determine the 
form of the general Killing tensor on E> n by considering it as an imbedded hypersurface 
in E n+1 . We also obtain algebraic conditions for the orthogonal integrability of the 
eigendirections of the CKTs. In Section 4 we formulate and solve the equivalence 
problems for CKTs defined on § 3 . Section 5 is devoted to applications to problems 
of classical mechanics that subsume the results of the previous sections. Section 6 
contains the conclusion. 
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2 Invariant theory of Killing tensors 



In what follows we formulate and solve the equivalence problem for the OSWs gener- 
ated by CKTs defined on S 3 within the framework of the invariant theory of Killing 
tensors (ITKTs), which is a natural extension of the classical invariant theory of 
homogeneous polynomials (see [H] for more details and relevant references). Thus, 
the study of OSWs is based on algebraic and geometric properties of Killing tensors 
two-tensors in which invariant theory comes into play as the natural link between al- 
gebra and geometry. Moreover, recently it has been shown explicitly (see [16] for more 
details and references) that in the case of Killing two-tensors generating OSWs the 
study can be naturally cast into the general setting of Cartan's geometry [5j EJ [12] . 

Let (Ai,g) be an n-dimensional (pseudo-)Riemannian manifold of constant cur- 
vature. 

Definition 2.1. A (contravariant) Killing tensor of valence p defined in (M.,g) is a 
symmetric (p, 0) tensor field satisfying the Killing tensor equation 

[K,g]=0, (2.1) 

where [ , ] denotes the Schouten bracket [32]. When p = 1, K is said to be a Killing 
vector field (infinitesimal isometry) and ( 12.1 ft reduces to 

C K g = 0, (2.2) 

where C denotes the the Lie derivative operator. 

Since the Schouten bracket [ , ] is R-bilinear, the set of solutions to the system of 
overdetermined partial differential equations (PDEs) given by (12. ip forms a vector 
space over M. Furthermore, since {Ai,g) is of constant curvature, the dimension of 
such a vector space is maximal (see the relevant references in [H] for more details). 
In what follows, we shall use the notation )C P (A4) to denote the vector space of 
valence p Killing tensor fields defined on M.. 

Remark 2.2. Note that the equation ( 12.11) can be equivalently rewritten in the more 
familiar (covariant) form as follows: 

% i P ;i J+1 ) = 0> (2.3) 

where K ix jip denotes the covariant components of the Killing tensor K, ; the co- 
variant derivative with respect to the Levi-Civita connection defined by g, and (. . .) 
symmetrization over the enclosed indices 
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As is well-known the Killing tensors defined on spaces of constant curvature are 
sums of symmetrized tensor products of Killing vectors forming finite-dimensional 
vector spaces. Denote by d the dimension of the vector space K, P (M.) and recall that 
d is given by 

Therefore the general element of )C P (A4) may be represented by d arbitrary param- 
eters a 1 , . . . , a d with respect to a given basis. Let G denote the isometry group of 
(A4,g). Our next observation is that the group G acting on M. induces, via the 
pushforward map, a linear transformation of 1C P (M.) which defines a representation 
of G [23]. Moreover, the action G O /C 2 (.M) is not transitive. These observations 
are the crux of the invariant theory of Killing tensors, allowing us to formulate the 
canonical forms and equivalence problems for the Killing tensors defined on spaces 
of constant curvature. 

Of particular importance for applications are the elements of K, 2 (Ai) enjoying 
additional geometric and algebraic properties. More specifically, let K G K?(M) be 
such that its eigenvalues are pointwise simple and real, and the eigenvector fields are 
normal (orthogonally integrable). Such a Killing tensor K is called a characteristic 
Killing tensor (CKT) [SJ[2]. The remarkable property of CKTs is that their eigen- 
values and eigenvectors generate OSWs, which are n foliations of the space M. that 
consist of (n — l)-dimensional hypersurfaces orthogonal to the eigenvectors of the 
CKT in question. 

The action G O JC 2 (A4), thus defined, foliates the vector space JC 2 (J\A) into the 
orbit space K?(A4)/G. These orbits represent Killing tensors of valence two that 
share the geometric (and algebraic) properties that are equivalent modulo the group 
action G O K?(M). Of particular importance in applications are the orbits that 
correspond to the CKTs belonging to the vector space K?(M). It must be noted at 
this point that, in general, the topology of the orbit space )C 2 (A4)/G (or {JC 2 (M) x 
M)/G) is far from being trivial which makes the problem of invariant classification of 
the orbits fairly complicated. The problem consists of two "subproblems," namely the 
canonical forms problem and the equivalence problem which can briefly be formulated 
in this setting as follows. 

1. Canonical forms problem: Consider the action G O 1C 2 {M). The problem 
is to determine the number of inequivalent orbits corresponding to the CKTs 
defined on (Ai,g) as well as the canonical forms representing each of them. 

2. Equivalence problem: Consider again the action G O K 2 (M). Let K G 
K 2 (M). First, the problem is to deterimine whether or not K is a CKT. If 
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the answer is "yes" , the main problem is to determine the corresponding orbit 
in the quotient space K?(Ai)/G that the Killing two-tensor in question K 
belongs to. Finally, we also want to determine the moving frames map pi [H)] 
that maps K to its respective canonical form. 

Recall that Eisenhart [8] outlined the solution to the canonical forms problem for 
CKTs formulated above and solved it for the case At = E 3 . In [16] (see also the 
relevant references therein) Horwood et al reformulated Eisenhart 's approach in the 
language of the Cartan geometry (which Eisenhart employed implicitly). Note that 
the canonical forms problems for the cases Ai = S 3 and At = H 3 were solved by 
Olevsky in [2H], while the case A4 = M 3 was treated by Horwood and McLenaghan in 
[To] . Moreover, these solutions have been used to solve the corresponding equivalent 
problems. Thus, Horwood et al [H] and Horwood [T7j employed two different meth- 
ods to solve the equivalence problem for the CKTs defined on E 3 , while Horwood at 
al [16] solved the equivalence problem for the case At = M 3 . 

In order to set the stage for our theory and in what follows solve the equivalence 
problem for )C 2 (S), we now employ the fundamental ideas from Cartan's approach to 
geometry. Observe first that Ai ~ G/H, where G is the isometry group of Ai and H 
is a closed subgroup of G. Thus, for example, S> 3 ~ SO(4)/SO(3). Next, we observe 
that in view of this identification the homogeneous space G/H can be treated as the 
base manifold in the tautological principal bundle projection tti : G — y G/H ~ Ai 
(i.e. G is the principal if -bundle over G/H). Similarly, consider the vector bundle 
lC p (At) x Ai, then 7r 2 o7r 2 : )C p (Ai) x At — > At ~ G/H is the vector bundle projection 
with the same base manifold. Upon noticing that the transitive action of the isometry 
group I(M) = G in Ai yields the non-transitive action of G in )C p (Ai), where G acts 
as an automorphism, we consider next the orbit space (K, p (Ai) x Ai)/G leading to 
the third projection 7r 3 : JC p (Ai) x At — > (K. p (Ai) x Ai)/G, having the structure of a 
principle G-bundle with JC p (Ai) x Ai as the total space. Following [IB] , we introduce 
the lift f : (K, 2 (Ai) x Ai)/G — > G, so that the following diagram commutes: 

G >G/H~M (2.5) 



/ 



7T 2 



{K P {M) x M)/G^ r lC p {M) x Ai 

Geometrically, the existence of / is equivalent to the existence of a cross-section 
through the orbits of the orbit-space (IC p (At) x Ai)/G. The intersections of such a 
cross-section with the orbits are the corresponding canonical forms, and their coordi- 
nates are covariants (invariants) which play an important role in the considerations 
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that follow. 

Consider now the case of p = 2. 

Let K G K, 2 (A4) be a CKT at a non-singular point x e M. (i.e. the eigenvalues 
of K are all real and distinct at x). Indeed, K gives rise to a quasi-ort honor mal 
frame E KjX (M) of eigenvectors {e 1; . . . , e n } of K at x G M., which is also a quasi- 
orthonormal basis for T X (M). Denoting E(M) as the corresponding bundle of frames 
generated by CKTs in A4, it follows that (E(Ai), A4, n 2 ) defines an (oriented) quasi- 
orthonormal frame bundle, where 7T2 : E(M) — > A4. The fibres 7r2~ 1 (x) correspond to 
sets of all possible quasi-orthonormal frames at x G M. generated by the eigenvectors 
of CKTs. Finally, this arrangement leads to the fibre bundle projection 7f 3 : )C 2 (A4) x 
M. — > E{M). Accordingly, our diagram (12. 5p now assumes the following form: 

G >G/H~M (2.6) 




(JC 2 (M) x M)/G K?(M) x M E{M) 

The existence of such a lift / is assured by the fact that G acts transitively on the 
bundle of frames E(M) for a given CKT K G JC 2 (Ai) and the classical Cartan 
lemma [12J: 

Lemma 2.3 (Cartan). Suppose that (p is a Q-valued one-form on a connected (or 
simply connected) manifold M.. Then there exists a C°° map F : M. — ^ G with 
F*u = if iff 

dip — ip A (p, 

where u is the Maurer- Cartan form on G. Moreover, the resulting map is unique up 
to left translation. 

Furthermore, we define the map F : E(Ai) — > G to be F = f o 7r 3 o tt^ 1 . Clearly, 
in this case G may be identified with the set of frames associated to the CKT K 
in question, or alternatively, cross- sect ions of the fibration G — > G/H over K. In 
view of Lemma 12.31 the Maurer-Cartan form on G can be restricted to this choice 
of frames and thus produce a complete set of invariants (covariants) to solve the 
equivalence problem, which one can now solve for the orbit space ()C 2 (A4) x M)/G 
(or fC 2 (M)/G) for CKTs using the classical calculus of differential forms. More 
specifically, the problem of invariant classification of the orbit (s) generated by a CKT 
K G K?(A4.) reduces to fixing a quasi-orthonormal frame of eigenvectors {ei, . . . , e n } 
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and considering in the frame the corresponding Cartan structure equations 

de a + oj\ A e b = T a , (2.7) 
du>\ + u> a c A u> c b = 0%, (2.8) 

together with the Killing tensor equations for the components K ab of K ( 12.31) and 
the integrability conditions 

e a Ade a = (no sum). (2.9) 

In these equations, uj a b = T cb a e c are the connection one-forms, T a = \T a bc e b Ae c are 
the torsion two-forms, @ a b = \R a hcd e c Ae d are the curvature two-forms, {e 1 , . . . , e n } 
is the dual basis of one-forms, the connection coefficients T cb a correspond to the Levi- 
Civita connection V (and hence T a = in (12. 7p ) and R a bcd are the components of the 
curvature tensor. Note that (12. 9p are the integrability conditions for the normality of 
the eigenvectors e a . We also note that, with respect to this frame, the components 
of the metric g and CKT K are given by 

0a6 = diag(ei,...,e n ), K ab = diag(eiAi, . . . , e n X n ), (2.10) 

respectively, where e a — ±1, a — 1, . . . , n, and A a , a = 1, . . . , n, are the eigenvalues of 
K . The differential invariants characterizing the orbits in question are determined 
by the connection one-forms u) a b that are found from a fixed quasi-orthonormal frame 
{ei, . . . , e n }. More specifically, solving the Killing tensor equation (12. 3 j) for K ab in 
this case modulo the integrability conditions ( 12. 9ft will produce a set of canonical 
forms (thus solving the canonical forms problem) corresponding to the orbits, while 
finding the connection one forms uj a b will provide a means to distinguish between 
the orbits, thus solving the equivalence problem (or, rather, most of it). However, 
at this point we would like to point out that solving the equivalence problem in this 
way (i.e., by thus finding a complete set of differential invariants) may be extremely 
challenging computationally. Instead, one can make use of the fact that the group G 
acts transitively on the bundle of frames and try to solve it "in the group" , employ- 
ing an algebraic approach to Cartan's method of moving frames j£l [TOj [29] . More 
specifically, via the tensor transformation laws, we can determine the action of G on 
the parameters ai, . . . , a<j (and the local coordinates of M. if necessary) that deter- 
mine the vector space K?(M) (or the product space K?(M) x M) and then find a 
set of complete invariants (covariants) as algebraic functions of the parameters (pa- 
rameters and local coordinates). This approach (analogous to the classical invariant 
theory of homogeneous polynomials) is also more preferable from the applications 
point of view, because the Killing two-tensors in Classical Mechanics (for example) 
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normally appear in terms of local coordinates on .M, rather than in the frame of its 
eigenvectors as in (12.1 Oft . 

The concept of a web symmetry, which signifies that the web is invariant under 
at least a one-parameter group of isometries, is another interesting geometric conse- 
quence stemming from the properties of covariants (eigenvalues) of CKTs. Suppose, 
for example, K is a CKT in K?(M) having n — 1 functionally independent eigenval- 
ues Xi, i — 1, . . . , n — 1. It follows that there exists a vector field V G K}(M) (V 
is necessarily a Killing vector), such that £v{K) = 0. Indeed, the integral curves of 
V are given by the common level sets 

n-l 

f^|{Aj = const}. 

i=i 

Now we can be more specific about solving the equivalence problem outlined 
above. Let K e )C 2 (Ai). We first verify whether or not the Killing two-tensor K 
in question is a CKT. This requires verifying that the eigenvalues of K are real and 
distinct, and that the eigenvector fields of K are normal. The latter step can be 
accomplished by verifying the vanishing of the Haantjes tensor (H condition) [13] 
(see [18J for more details): 

H k (X, Y) := KN k (X, Y)+N k (KX, KY)-K(N k (X, KY)+N k (KX, Y)) = 0, 

(2.11) 

where the (1, l)-tensor K = Kg 1 , N k is the Nijenhuis tensor [27] and X, Y are 
arbitrary vector fields on A4. Thus, the vanishing of the Haantjes tensor defined by 
(12. lip is equivalent to the Killing tensor K being a CKT. Alternatively, the formula 
(12. lip can be given in index form as follows: 

W jk = N\ m K l 3 K m k + 2N e m[j K m k] K i £ + N l ' jk K m t K k m = 0, (2.12) 

where K l j denotes the components of the (1, l)-tensor K and [. . .] denotes skew 
symmetrization over the enclosed indices. 

One can also employ the Tonolo-Schouten-Nijenhuis (TSN) conditions to verify 
the normality of the eigenvectors of a symmetric tensor field K of valence (0,2). 
Indeed, K with real distinct eigenvalues has normal eigenvectors iff the following 
conditions are satisfied: 

N\ijgk]e = 0, 

N%jK k]£ = 0, (2.13) 
N [ijK k ] m K m e = 0, 
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where N l jk are the components of the Nijenhuis tensor [27] of K l j defined by 

N' jk = K\K e m + K%K\ U . (2.14) 

In the index-free form, formula ( 12 . 141) can be obtained from the formula (12.111) by 
replacing the operator N^(-, ■) with the Lie bracket [ , ]. 

Remark 2.4. If the eigenvalues of K are real and distinct, the H condition and 
the TSN conditions are equivalent. If the multiplicity of some the eigenvalues of K 
is greater than one, the TSN conditions no longer apply. However, in this case the 
following more general result [13] needed in Section 4 is true: the distributions defined 
by the eigenspaces of K in a Riemannian space are orthogonally integrable if and only 
if the H condition is satisfied. Recall that orthogonal integrability signifies that the 
distribution defined by the orthogonal complement of the eigenspace is integrable. 

If K E K?(M) is a CKT, our next step is to determine which orbit in )C 2 (A4)/G 
it belongs to, noting that the corresponding canonical form for the orbit in question 
is known from solving the canonical forms problem. This problem can in principle 
be solved by employing a number of mathematical tricks and tools, although it is 
hard in general due to the fact that the group action G O /C 2 (.M) is not regular and 
thus the orbits have different dimensions. 

To solve the problem, one can employ algebraic invariants of the group action 
G K?(M). These are the invariants of the group G acting in the parameter space 
of the vector space /C 2 (.M), or more specifically, functions of the parameters that re- 
main unchanged under the induced action of the group G (see, for example, [HI [32] 
for more details). It must be noted however, that the problem can normally be 
solved in relatively simple cases using only algebraic invariants, that is when one 
has to deal with few types of orbits (e.g., dim.M is small). Recall that algebraic 
invariants of the CKTs were first introduced by Winernitz and Fris [32] and later 
rediscovered in McLenaghan et al [22] in a different context. More generally, one has 
to recover more information about the orbits and their degeneracies. One way to deal 
with the problem is to "improve" the group action by considering it in the extended 
space K?(M) x M.. After appropriate computations, this yields algebraic covariants 
of Killing tensors, which are functions of both the parameters of the vector space 
/C 2 (A4) and coordinates of Ai. The covariants of Killing tensors were introduced 
by Smirnov and Yue in [33J and have been successfully employed to solve a number 
of equivalence problems (see, for example, Horwood [T7]). The degeneracies of the 
orbits in K, 2 (Ai)/G (or ()C 2 (Ai) x Ai)/G) manifest themselves in the orbits having 
different dimensions - a fact which for our problem has appropriate group theoretical, 
geometric and algebraic interpretations. Thus recall that in general the dimension of 
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an orbit is determined by the formula dim G = dim X + dim G x , where G is a group 
acting on M, x e M, O x is the orbit through x, and G x is the isotropy subgroup of 
G through x. The existence of isotropy subgroups indicates the degeneracies of the 
orbits (i.e., their dimensions drop), which in our case of Ai = JC 2 (M)/G from the ge- 
ometric perspective is equivalent to the CKTs corresponding to the degenerate orbits 
admitting web symmetries (see below). Algebraically this means that the CKTs cor- 
responding to the degenerate orbits have functionally dependent eigenvalues. Note 
that if all the eigenvalues of a CKT K are functionally independent they can be 
taken as the orthogonal coordinates (see [16] for more details). Degeneracies of the 
orbits can also be characterized by singular points of the CKTs in question, that 
is the points where the eigenvalues coincide. Unfortunately, most of the orthogonal 
coordinate webs used in applications come from the CKTs with degeneracies. 

Solving the equivalence and canonical forms problems described above forms the 
mathematical foundation for applications arising in Classical Mechanics, in particular 
in the Hamilton- Jacobi theory of orthogonal separation of variables for the natural 
Hamiltonians defined in spaces of constant curvature (see [16] for relevant references 
and details). Let a Hamiltonian system be defined by a natural Hamiltonian of the 
form: 

H(q,p) = \g lj (q)p lP , + V(q). (2.15) 

The following theorem (Benenti [2j), which is a generalization of the corresponding re- 
sult due to Eisenhart concerning geodesic Hamiltonians [8], establishes an important 
link between orthogonal separation of variables in the associated Hamilton- Jacobi 
equation for the Hamiltonian system defined by (12.151) . 

Theorem 2.5. A Hamiltonian system defined by ( TJS5P is orthogonal separable iff 
there exists a CKT K such that 

d(KdV) = 0, (2.16) 

where the (1, l)-tensor K = Kg 

In light of Theorem 12.51 the canonical and equivalence problems formulated above 
can be reformulated in the language of orthogonal separation of variables for natural 
Hamiltonian systems defined by (I2.15P in spaces of constant cutvature as follows: 

(i) How many "inequivalent" coordinate systems afford orthogonal separation of 
variables in the corresponding HJ equation? 

(ii) If the answer to (i) is non-zero, how can one characterize intrinsically the 
coordinate systems that afford separation of variables in the HJ equation? 
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(iii) What are the canonical coordinate transformations 

{q\q\...,q n )^{u\u\...,u n ) 

from the given position coordinates of (12.151) to the coordinate systems that 
afford orthogonal separation of variables of the HJ equation? 

In what follows we employ some of the tools and techniques described above to 
solve the equivalence problem for the CKTs defined on S 3 , employing the results ob- 
tained by Olevsky [28] and Kalnins et al [20] who solved the corresponding canonical 
forms problem. 



3 Killing two-tensors on E> n C E ra+1 

To fully explore the algebraic and geometric properties of Killing two-tensors defined 
in a sphere, we employ the (Cartesian) coordinates of the corresponding ambient 
Euclidean space. In this view the vector space of Killing two-tensors /C 2 (§ n ) is viewed 
as a subspace of /C 2 (E n+1 ) whose elements are defined in terms of the Cartesian 
coordinates of E n+1 . 

Indeed, since every Killing tensor defined in a space of constant curvature is 
expressible as a sum of symmetrized products of Killing vectors, let us begin by 
defining basic Killing vectors of a Euclidean space E m in terms of the corresponding 
Cartesian coordinates. For our purposes, we also need to define the dilatational 
vector field D (Euler vector field). Thus, the translational Xi and rotational R^j 
Killing vectors and the vector D in terms of a system of Cartesian coordinates x % 
can be defined as follows: 

Xi = £j, R ij = 28%g tm x m X k , D = x l X u (3.1) 

where g^ = diag(l, . . . , 1), denotes the Euclidean metric of the Euclidean space in 
question and 6^ = SHS^, the generalized Kronecker delta. Note that the translational 
and rotational Killing vectors defined above form a basis for the Lie algebra /C 1 (E n ) 
of the group of rigid motions of E m . Next, the commutation relations among these 
vectors are given by 

[Xi, Xj] = 0, [Xi, Rjk] = 25j™g mi X£, [Xi, D] = X h ^ 
[Rij, Rke] = 4:S^ n 5l r e g mp R nr , [D, Rij] = 0. 
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The Killing vectors also satisfy the following algebraic identities (syzgies): 

X[i R jk] = 0, R i{j R kl] = 0, (3.3) 

where © denotes the symmetric tensor product. The general Killing vector thus has 
the form 

K = A i X i + B ij R ij , (3.4) 

where the constants A 1 , and B % \ called the Killing vector parameters, satisfy the 
symmetry relation 

BM = 0. (3.5) 

From (13. ip and (13. 4ft we obtain the covariant components of K with respect to the 
natural basis dx l : 

K i = A i + 2B ij x j . (3.6) 
Similarly, the general valence two Killing tensor thus has the form 

K = A ij X { © X 3 + B^Xi R jk + (jWRij R M , (3.7) 

where the constants A l \ B l ^ k , and C l ° : , called the Killing tensor parameters satisfy 
the following symmetry relations 

= 0, 

B i(Jk) = Q) B m = Qj (3.8) 

^j(ij)k£ g ^jijkl ^jklij ^ji[jk£] q 

A key observation is that the parameter set has the same symmetries as the 

Riemann curvature tensor. For this reason it is sometimes called an algebraic curva- 
ture tensor. From (13 .ip and (13. 7p we find that the covariant components of K with 
respect to the natural basis are given by 

K i:j = Aij + 2B {ij)k x h + 4:C ikje x h x e . (3.9) 

This result is consistent with that obtained in [23] obtained by the use of represen- 
tation theory. 

Now we may determine the form of the general Killing tensor on § n by considering 
it as an imbedded hypersurface in E m with m = n + 1 defined implicitly by the 
equation 

x-x = l, (3.10) 
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where x = (x 1 , . . . ,x n+1 ), and • denotes the Euclidean inner product. Let 

x = f(u\...,u n ), (3.11) 

be a local parametrization of S n . Then from (I3.10p and (13. lip we obtain by differ- 
entiation 

x ■ x a = 0, x af} ■ x = -g ap , (3.12) 
where x a = J^, (a^ = and 

9a/3 = x a ■ Xp, (3.13) 

denotes the pullback of the Euclidean metric to S n . The characterization is given by 
the following proposition: 

Proposition 3.1. The general Killing vector and general valence two Killing tensor 
fields on S n are given by 

K a = x^Kt, K a(3 = x^Kij, (3.14) 

where 

Ki = 2B ijX j , Kij = AC im x k x l . (3.15) 

Proof. The pullback to S n of the covariant derivative of any smooth tensor field 
Ki lt ,,, t i v defined on E n+1 is given by 

X a\i ■ ■ ■ J X a p X a p +i^-h,—,ip,ip+l ~ ^ a.i,...,a p ;ap+\ + P x ^!(a2,->«p-i3a P )a p +r (3.16) 

If K^...^ is any solution of the KT equation fl 2 . 3 j) . it follows that 

K(a 1 ,...,a p ;a p+1 ) + PX t, Ki(a 1 ,...,ap- 1 9a p a p+1 ) = 0. (3-17) 

We conclude that K will define a Killing tensor field on § n if and only if 

£ l -^i(ai,...,ap-i3'Q !p a p+ i) = 0. (3.18) 

This condition is clearly satisfied if 

■r J K ji: : = 0, (3.19) 

on E n+1 . For the cases p = 1, 2 this condition takes the form 

x*Ki = 0, x j K {j = 0. (3.20) 
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It follows from (j33J, (l3l)]) . (ETSjl . and that 



A i = 0, ^ = 0, 5 fo>=0 . (3.21) 

These conditions imply that and ifjj have form given by ( 13.151) . We note that, 
with the use of ( 13. 2ft . the conditions ( 13 . 2 1 f) may be written invariantly as 

[D,K] = 0. (3.22) 

It remains to demonstrate that (13.141) and (13.151) indeed define the general Killing 
vector and general valence two Killing tensor on S n . This may be achieved by a 
dimensional argument. From (12 .4p it follows respectively for p — 1, 2 that 

d=-n(n + l), -^-n(n + l) 2 (n + 2). (3.23) 

For the case p = 1, the number d is identical to the number of independent compo- 
nents of the Killing vector parameter Bij, while for the case p = 2, d is identical to 
the number of independent components of the algebraic curvature tensor CV,h in a 
(n + l)-dimensional space. □ 

Remark 3.2. Proposition ^. H is closely related to a result proved by Delong [7] which 
states that the valence p Killing tensors (considered as functions on the cotangent 
bundle of E n+1 ) that are in involution with the Euclidean distance function and the 
function x l pt, where pi denotes the canonical momenta, are Killing tensors on S n . 
Delong's conditions are equivalent to (13.201) and (13.221) respectively, only one of which 
seems to required to obtain the result. Furthermore, our proposition goes beyond his 
by giving the explicit form of the general valence one and two KTs on S n in terms 
of the appropriate KTs in the ambient space E n+1 . Our proposition may be easily 
generalized to valence p Killing tensors by the use of the results of Horwood |17j . 

Remark 3.3. Note that the image of the contravariant metric metric g al3 on S n with 
respect to spherical coordinates under the pushforward induced by the parametriza- 
tion gives rise to a Killing tensor C called the Casimir tensor with the property that 
it commutes with every element of the vector subspace of /C 2 (E" +1 ) given by (13.151) . 
In view of the above C G /C 2 (E™ +1 ) is degenerate as an element of the vector subspace 
determined by 4Cikjex k x e . We also note that C is a CKT for the spherical web on 
S n . The explicit form of the Casimir tensor is given by 

C = kg^Ra Rkt, (3.24) 
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or with respect to the natural basis by 

Cij = kg i[:j g e]k x k x e (3.25) 

where k is some constant. It may be shown by direct calculation that [C, K] = 0, 
which verifies the above mentioned commutation property. 

The next step is to determine the KTs among those given by f)3.14p and (" 13 . 1 5[) 
that are CKTs, that is which have point- wise distinct eigenvalues and normal eigendi- 
rections. To effect this determination we use the fact the eigenvalues of K a p are also 
eigenvalues of and that the pushforward of the eigenvectors of K a p are eigenvec- 
tors of K^. These results follow by contracting both sides of the eigenvalue equation 

KyX* = \ gij X\ (3.26) 

with x l a , which yields 

K aP X? = Xg aP X^ (3.27) 
where X- 7 in (I3.26f) is the pushforward of X@ 

X j = x ] p X p . (3.28) 

We conclude that if has distinct eigenvalues then so does K a p. The converse 
is not true since if K a p has a zero eigenvector then by (I3.20p will have zero 
as a repeated eigenvalue. However, the remaining eigenvalues will be non-zero and 
distinct. We also observe by A3. 12[) that the eigenvector x l corresponding the the 
zero eigenvalue of pulls back to the zero vector on § n . 

We next show that the pullback of a normal co-vector field on E n+1 is a normal 
vector field on S n . To see this consider any co- vector field E^ on E n+1 which satisfies 
the normality condition 

E [id E k] = 0. (3.29) 

We note that (I3.29p is the component form of the integrability condition f 12 . 9 j) . An 
easy computation shows that 

EijEkxi^x* = E a> pE 7 - x^EiEy. (3.30) 

It follows from fl3T29|) and fl3T30|) that 

E[a,pE 7 ] = 0, (3.31) 

which implies that the pullback E a is a normal vector field. The above results imply 
that the pullback of a characteristic KT on E n+1 is a characteristic KT on § n . In 
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view of the above results we are able to study CKTs on § n as pullbacks of CKTs on 
E n+1 . 

We observe that the use of ( \2.9\\ to determine whether the eigenvector fields 
of K are normal, is impractical since it requires the explicit determination of the 
eigenvector fields, an intractable problem. However, alternatively we may utilize 
either the TSN conditions ( 12.131) or the H condition defined by ( 12.121) . We first 
consider the TSN conditions for a KT of the form ( 13.151) . These conditions impose 
the corresponding algebraic conditions on the algebraic curvature tensor C: 

C e ( pq [iCjk]r)e = 0, (3.32) 

C£( pq m C i r [ijCf : ]st)m ~ 2C^(p 9 [j(7j \ r \ m Ck]st)m = 0, (3.33) 
°W(pg ^ r\n\s^ t[ij { ^'k]uv)m > ^i{pq ^^rs^ j\t\ ^k]uv)m ' 

%C^ pq m C\ ri \ TS ] i iCj n \ t \Ck}uv)m — 0, (3.34) 

where | . . . | denotes exclusion of the enclosed indices from the symmetrization pro- 
cess. The first TSN condition (13.321) implies 

C£( pq rn C l ' r [ijC k]st)m + 2C 'e(pq[iC 'j '\ r \ m C 'k]st)m = 0- (3.35) 

This equation together with the second TSN condition (I3.33j) yields the two equations 

Ce( pq m C e r [ijCk] s t)m — 0, (3.36) 

Ct(p q [iC/\ r \ m Ck]st)rn — 0. (3.37) 

These equations replace the second TSN condition (I3.33[) . We now show that the 
third TSN condition is a consequence of the first and second conditions. The cyclical 
identity implies 

^l(pq ^^rs^j \t\ ^kjuv^ — ~ { ~'£(pq <— ' |n|rs [i 1 -^ j |i | t> k]uv)m~T ^ i(pq ( -'|n|rs[i ( - y j \t\ ^k\uv)m- 

(3.38) 

From the above equation and (I3.37P we obtain 

Cl(pq m C\n\rs[iCj\t\ ni Ck]uv)m = 0, (3.39) 

and 

e(pq ^\n\rs[i^j \t\ ^k]uv)m — ^l{pq ^\n\rs[i^j \t\ { ~ y k]uv)m- {O.W ) 
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Now from the first TSN condition and (13.401) we find 

Crrifil fin fi or 1 in ft fi n £ ft _ r\ (o A'\\ 

l{pq ^ r\n\s^ t[ij^k]uv)m ~ ^ l(pq ^\n\rs[i^j \t\ ^fc]™)?™ — U. I'J-^-'-J 

On the other hand (I3.36j) implies 

fi nifit fin fi I nf m fi fi n £f _ r\ (o /|Q\ 

^£(pq rlnls*- 7 t[ij t> k]uv)m T ZLs£( p q Ls\ n \ rs [iLsj \t\ ^k]uv)m — U. [O.^iZ) 

These equations together imply 

C '(.(pq" 1 C ' r \ n \ s C n 't[ijC 'k] U v)m = 0, (3.43) 

Cz(p q m C\ ri \ rs [iCj n \ t \Ck}uv)m = 0. (3.44) 

Finally (ESSj), (EE]), and fjTO^It imply that the third TNS condition (EgJ) is iden- 
tically satisfied. 

Remark 3.4. This result was first proven in E 3 by Czapor [H] using Grobner basis 
theory and computer algebra. Schobel [31] extended the result to n-dimensional 
spaces of non-zero constant curvature using representation theory of the symmetric 
group. However, it seems that our proof, based on standard indicial tensor algebra, 
is simpler and more concise. 

Substituting the Killing tensor ( 13. 151) into the H condition (12 . 121) we obtain the 
following condition on the coefficients: 

Afi kfim fi nfii I r\/~i kfin fi mfil 

^^(.{pq ^ rs\i^ j\tu t — ' v)mn > '£(p|m| ^ qr\v^ j]st ^ uv)n 

zfi kfim fi nfii , fi kfim fi i nfi 

d^ilpq ^ rs[j t ~ / j]|m|i ^ uv)n "T ^ l(pq ( - / rs[i^j] t ^\n\uv)m 

I fi kfim fi nfi I Q/^ kfim fin fi£ 

~r'~ y £(pq rs[i^j]tu *~^\n\v)m ^^l.(pq *~ y r[ij'|^ stlml*- 7 uv)n 

— 1Cl{pq C rs[iCj]t\ m \ C uv)n = 0. (3.45) 

Using this condition we can verify whether or not a given Killing tensor defined on 
§ n is characteristic. It may be shown that that the first and second TSN conditions 
(I3.32p and (I3.33j) imply that the H condition (I3.45j) is satisfied. 
We next examine the pullback of the quadratic first integral 

K = K ij p iPj + U. (3.46) 

The vanishing of the Poisson bracket implies that Kij satisfies (12. 3 p for p = 2, and 
that 
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where V is the potential function. The pullback to S n is given by 

dx { 3 13 1 dx l ' 



£4 = 9^4x^x1 (3.48) 



which may be written as 

We conclude this section by studying the action of SE(m) on the Killing vector 
and Killing tensor parameters [17J. The action SE(m) O E m is given by 

x l = A l jX j + 5\ (3.50) 

where A l j G SO(m), 5 l G E m , and x % denote the transformed Cartesian coordinates. 
This tranformation induces by ( 13. ip the following transformation of the Killing vec- 
tors: 

Xi = A J \Xj, Rij = h k i A l jR kt + fXij k X k , (3.51) 

where 

Pi/ = 25^g mn A k t 5 m . (3.52) 
It follows that the Killing vector and parameters transform as 

A* = A^A*' + fi^B^, B ij = A i k A j t B M , (3.53) 

while the valence two Killing tensor parameters transform as 

I« = A\A^A kt + 2A^ m ^B klm + iikivmn j C klmn , 
B^ k = A\AJ m A k n B £mn + ^ mn l A\A k q C mn ^, 

C^ kt = A i m A j n A k p A t q C mnpq . (3.54) 

In order to make our formulas more compact we introduce a multi-index notation 
where any upper-case index represents a pair of skew-symmetric lower-case indices. 
Thus Ri will represent Rij. Using this notation we may rewrite (I3.5ip as 

Rj = A J I R J + ^ I k X k , (3.55) 

where A 1 k = AVA^i, represents the second compound of A l j. With the same 
notation (I3.53P reads 

A 1 = A l jA j + fij l B J , B 1 = A 1 jB J , (3.56) 
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whereas ( 13.541) reads 

A ij = A\A J e A ke 
B iJ 



2A (l 



A\A J M B m 



iKL 



k fi L j) B kL + ^^C 1 
md + fi L l A J mC lm , 
C IJ = A 1 L A J M C LM . 



(3.57) 



We are now in a position to give an invariant classification of the translational 
Killing vectors on E m under the action of SE(m). The invariant needed to effect this 
classification is T\ = B^B 1 ^ . 

A translational Killing vector in E m is invariantly defined by the condition T\ = 0, 
which implies B l i = 0. Thus by (13. 4 p a translational KV has the form 



K = A l X, 



(3.58) 



where it is assumed that not all the A 1 are zero so that K is non-trivial. It follows 
from the transformation formula (I3.53P that, by an appropriate choice of A l j e SO(m) 
we can transform (13.581) to the form 



K 



(3.59) 



for some A 1 ^ 0. 

The case when X\ ^ 0, characterizes the non-translational Killing vectors the 
analysis of which will be left to a future paper. 

However, we are in a position to classify the rotational Killing vectors in E 4 
given by (13.151) that define KVs on S 3 . For this we need the additional invariant 
X 2 = €ijkiB 13 B kl . We utilize the property [11] that the skew-symmetric matrix B l i , 
assumed to be not the matrix, may be transformed to the following canonical form 
by an element of 0(4): 



B' 



( 










o\ 




-«1 
























\ 








-«2 


0/ 



(3.60) 



where K\, k% G R. If X\ ^ and X 2 = 0, then k\ + k\ ^ and = 0. Without 
loss of generality we may set k 2 = 0, which implies that the KV has the form 



K = b 12 Ru, 



(3.61) 



where 612 = 2/«i is some constant. We now assume X 2 7^ 0, which implies that 
7^ 0. In this case the matrix B l i is rank four from which it follows that the KV 
has the form 

K = b 12 R 12 + b u R u , (3.62) 
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where &12 = 2«i and b±2 = 2k>2 are some constants. The equations ( 13.611) and ( 13.621) 
via Proposition 13.11 give the two possible canonical forms for the rotations on S 3 . 

4 The equivalence problem for Killing tensors on 

To solve the equivalence problem outlined in the preceding section for the CKTs 
on § 3 we employ as in the Cartesian coordinates of the ambient Euclidean space 
E 4 D § 3 . In terms of these coordinates, the general Killing tensor of S> 3 is given by 

K = 4:C ijke Rij Ru, (4.1) 

where 

Rij = 8 m ijkX X m , 5 m ijk = &T9jk — 8J l gik, (4.2) 

i,j,k,£,m = 1,...,4, 

x 1 , . . . ,x A are Cartesian coordinates of E 4 , Xj = 9$, <7y are the components of the 
metric and is the symmetric tensor product. Note that ( 14.21) defines the six gen- 
erators of so(4), the Lie algebra of the isometry group 5*0(4) of § 3 (no reflections). 
Therefore the RHS of (14. ip represents a vector space of Killing tensors on § 3 . Fur- 
thermore, the dimension d of this space is determined by the tensor C^ M which has 
the same symmetries as that of the Riemann curvature tensor, hence d = 20. 

The special orthogonal group SO (4) is a Lie subgroup of the orthogonal group 
0(4), consisting of all orthogonal matrices A with positive unit determinant. The 
transitive action of SO (4) on E 4 can be specified by 

x l = A^x 4 , 

where A^- G SO (4) and x % denote Cartesian coordinates. This, in turn, induces the 
following transformation 

Ru = A\A £ jRki, & = A\A\K M (4.3) 

on the Killing vectors (14. 2p and Killing tensors ( 14. ip of E 4 . At the same time, this 
action induces the following transformations 

B ij = A\A\B H , C ijke = A 1 p A j q A k r A £ s C pqTS (4.4) 

on the Killing vector and Killing tensor parameters. 
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To make our formulas more compact, we will once again adopt a multi- index 
notation (i.e., Ri will represent fly) . Using this notation we may rewrite the first 
equation of ( 14. 3ft as 

Ri = \ j iRj, (4.5) 

where X 1 K = A'^A-^, represents the second compound of X l j. With the same notation 
(1441) reads 

B 1 = A J j5 J , C* /J = A / x A J lC xl . (4.6) 

To obtain the invariants and covariants of the group action SO (4) O /C 2 (§ 3 ), we 
apply the theory developed in [17] . In particular, by taking contractions of products 
of the general Killing tensor, the Euclidean metric, and the coefficient tensor, we can 
obtain a complete set of SO(4)-covariants and invariants. Using matrices K and C, 
where (K) l j = K 1 ^ and (O) 1 j = C lJ ki , as well as the trace operator "Tr", we can 
express these covariants Ci and invariants I, as follows 

Ci = Tr(K), C 2 = Tr(K 2 ), C 3 = Tr(K 3 ), C 4 = 9ij x l x\ 

I 1 = Tr(Q, l 2 = Tr(C 2 ), X 3 = Tr(C 8 ), ... ,J 14 = Tr(C 14 ), 

where, for example, Tr(K) = K l i and Tr(C) = C 1 1 . 

As we stated in Section [2, the action SO(4) O /C 2 (§ 3 ) foliates the vector space 
/C 2 (§ 3 ) into the orbit space /C 2 (S 3 )/SO(4). Each orbit is represented by a canonical 
form, and the solution to the equivalence problem requires determining such canon- 
ical forms as well as a classification scheme for finding which orbit a given CKT 
belongs to. As we mentioned in the introduction, the canonical forms problem has 
been solved. Please refer to the appendix for a list of six canonical forms for the 
orbit space /C 2 (§ 3 )/SO(4). 

Let us now develop a classification scheme for the orbit space /C 2 (S> 3 )/SO(4). The 
set of invariants and covariants listed above could be used to try and classify the orbits 
of these CKTs. This approach was successfully implemented in the solution to the 
equivalence problem of Killing tensors defined on E 3 [17], although the calculations 
were quite cumbersome. Indeed, the difficulty in this approach lies in finding certain 
linear combinations of the above invariants and covariants which distinguish between 
the orbits. A different approach which has proven to be more efficient and successful 
for the problems with small numbers of orbits is the method of web symmetries. 
The central idea of this method is to use the symmetry properties of the associated 
orthogonal separable web of a canonical CKT to characterize its orbit. In Section 2 
we stated that the web symmetries of a Killing tensor K are generated by Killing 
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vectors on the manifold. Thus, to determine all of the symmetry generators of a 
given web we impose the following condition 

C V K =0 (4.7) 

on the Killing tensor K defining the web, using the general Killing vector V of our 
manifold. 

In what follows we demonstrate the surprising result that the six CKTs of S 3 can 
be classified based purely on the symmetry properties of their associated webs. To 
this end, we are interested in obtaining the symmetry properties of a web before we 
impose the spherical constraint, which will yield additional web symmetries for the 
CKTs. Visually, this corresponds to capturing all of the symmetry properties of a 
web before it is intersected with the surface of S> 3 . 

To achieve this, we impose condition ( 14. 7p on each of the six CKTs using the 
general Killing vector of the ambient space E 4 . This will enable us to determine if a 
web is rotationally and/or translationally symmetry before it intersects the surface of 
§ 3 . Upon applying this method, we find that four of the six webs admit at least one 
rotational web symmetry. We can go even further by noting the number of rotational 
symmetries a CKT admits, which effectively divides the six canonical forms into three 
categories. Lastly, we find that two of the six webs admit translational symmetry, 
which provides the final distinguishing feature between each of the six webs. Please 
refer to Table H] for a summary of these results. 

Remark 4.1. In an application problem, it is possible that a given CKT K may 
have the Casimir tensor present. Specifically, 

K=aC + K- U 

where C is the Casimir tensor (I3.24p and K\ is a CKT. If Ki is translationally 
symmetric, then the addition of the rotationally symmetric Casimir tensor destroys 
this translational symmetry. Thus in order to determine all of the symmetries of K 
with or without the presence of the Casimir tensor, it is necessary to check the more 
general condition 

C v (K+aC) = 0, 

for arbitrary a. 

It is necessary to prove that the aforementioned symmetry properties of a Killing 
tensor are invariant under the action of SO (4). To do so, it suffices to solve the 
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Table 1: Symmetry classification of the six webs of S 3 under the action of SO (4) 



Category 


Symmetry 


Separable webs 


Generators 


I. 


2 rotations 


cylindrical 


Rl2 ? R-34 


II. 


1 translation & 1 rotation 


spherical 


R12, X 4 


III. 


1 translation 


spheroelliptic 


x 4 


IV. 


1 rotation 


elliptic-cylindrical I 


R42 






elliptic-cylindrical II 




V. 


none 


ellipsoidal 





Table 2: Invariant classification of Killing vectors on E 4 under the action of SO (4) 



Category 


Canonical form 




I 


R12 




II 


x x 






equivalence problem of Killing vectors of /C 1 (E 4 ) under the action of the group SO (4). 
To begin, we note that the general Killing vector of E 4 is given by 

K=A i X i + B I R I , 

where A 1 and B 1 denote the Killing vector parameters, and Xj and Rj are the Killing 
vector fields defined previously. The action of SO (4) on K induces the following 
transformations 

A = A\A\ B 1 = A I K B K , 
on the Killing vector parameters. Therefore, it follows that 

I 1 = B T Bj, 1 2 = A% 

are invariants in the orbit space /C 1 (S 3 )/SO(4). Using either of these two invariants 
it is possible to distinguish between two different types of symmetry generators. 
Please refer to Table H] for a summary of these results. We can conclude that the 
translational and rotational web symmetries as defined by the Killing vectors of 
/C 1 (E 4 ) are inequivalent under the action of SO (4). 

In addition to a classification scheme, a solution to the equivalence problem also 
requires a method for determining the moving frames map which identifies the group 
action required to return a given CKT to the canonical form of its orbit. On the 
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two-dimensional manifolds E 2 ,M 2 and S> 2 , algebraic formulas have been derived [221 
1241 IT~8] for determining the moving frame map of a given CKT. On E 3 and M 3 , a 
combination of web symmetry and eigenvalues and eigenvectors of the parameter 
matrices has be used to determine such a map [T4"l ITS] . In our case, however, the 
situation is complicated by the fact that our coefficient tensor, C ij , has order six 
when regarded as a matrix. As such, we will need to devise a different strategy for 
determining the moving frame map of a CKT on S> 3 . 

It has been noted in Section 3 from (13.81) that the coefficient tensor C 1 ^ has 
the same symmetries as the curvature tensor, and thus can be called an algebraic 
curvature tensor. In light of this property, let us lower the last three indices of C 1 ^ 
and contract on the first and third indices 

Ric = C l j^ = Rjg 

to obtain an algebraic Ricci tensor. The coefficient tensor for each of the six canonical 
forms listed in the Appendix can be contracted to define a canonical Ricci tensor in 
each case. The following proposition demonstrates that the Ricci tensor can be used 
to define the moving frame map for a given CKT. 

Proposition 4.2. A Killing tensor fl^. 1\ ) is in canonical form if and only if its Ricci 
tensor is in canonical form. 

Proof. Since the canonical form of the Ricci tensor is defined by the canonical form 
of K, the first direction is trivial. For the other direction, we prove by contradiction. 
Suppose the Ricci tensor of a Killing tensor K is in canonical form, but K is not. 
Since SO(4) acts transitively on the orbits of /C 2 (S 3 )/S'0(4), we can find a group 
action A G SO (4) which sends K to its canonical form K. In particular, the compo- 
nents of K transform according to (14. 3 P which induces the following transformation 

u jki ~ IV m 1Y j 1V k IV npq 

on the coefficient tensor C l j ke . At the same time, this action on C induces the 
following transformation 

Rj£ = A-^jA^iRmn 

on its Ricci tensor 1Z. Since the Ricci tensor of a canonical Killing tensor is necessarily 
canonical, we must have 1Z = TZ. This is a contradiction. □ 

According to Proposition I4.2| the moving frame map of a CKT can be constructed 
by determining the moving frame map of the corresponding Ricci tensor. Note that 
each canonical Ricci tensor can be represented by a diagonal matrix of order four. 
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Therefore, the determination of the moving frame map for the Ricci tensor is an 
eigenvalue-eigenvector problem for matrices of order four. Before we illustrate this 
technique with the application in the next section, we summarize our results in the 
following algorithm. 

1. Begin by substituting the potential into the compatibility condition (12.161) to 
determine the most general Killing tensor compatible with the potential. Using 
this Killing tensor, determine the subspace of CKTs. 

2. Next, we classify a CKT K by determining whether it admits any symmetry. 
Namely, impose the constraint 

C v (K+aCf) = 0, 

where a is an arbitrary parameter, C is the Casimir tensor, and V is the general 
Killing vector of E 4 , 

K= A l Xi + S / J R / . 

If K does admit symmetry, determine which type and the number of generators 
for each type. Consult Table H] to classify the CKT. 

3. To determine the moving frame map for K, find the Ricci tensor 1Z of the 
coefficient tensor. Diagonalize 1Z by solving the corresponding eigenvalue- 
eigenvector problem. The matrix A, which diagonalizes 71 defines the moving 
frame map. 

4. Finally, define the orthogonally separable set of coordinates corresponding to 
K by substituting A found in the previous step into the equation 

x* = A i j T j {u k ), 

where x % = T^{u k ) denote the canonical orthogonally separable coordinates 
corresponding to K. 

5 Application 

Consider the following natural Hamiltonian 

H = C i3 ViVj + - 1 

[x-yy 
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defined on S 3 , where C 8 - 7 denotes the Casimir tensor and x, y, z, w are cartesian coor- 
dinates of the ambient space E 4 . Using this Hamiltonian, we will now demonstrate 
how to apply the theory outlined in this paper. 

First, we impose the compatibility condition (I2.16P to obtain a family of Killing 
tensors which are compatible with the potential. Of this family, the following re- 
strictions on the parameters yields a subfamily of Killing tensors which satisfies the 
Haantjes condition ( 12. lip and generally admits 3 distinct eigenvalues: 



C1212 
C1323 



C3434, C1313 



C*2323, Cl414 



a 



c 



1313 



Cl212, Cl424 



c 



1212 



2424, 
Cl414 



Therefore we conclude that K must characterize at least one of the six orthogonal 
separable webs of S 3 . After a direct calculation, we find that K admits the following 
family of rotational Killing vectors 

3 3 3 3 

V = (C 3 Z + C&W)— + (CqW - C 3 Z )g^ + ( C 3V ~ C 3 X )g^ ~ ( C 6V + 

for arbitrary constants C3 and c$. Using the classification scheme outlined in Table 
H]we conclude that K characterizes a non-canonical cylindrical web. 

In order to determine the orthogonally separable coordinates for this Killing ten- 
sor, we need to determine the transformation which maps K to its canonical form. 
As discussed in Section HJ such a map can be constructed by diagonalizing the Ricci 
tensor of the coefficient tensor. Contracting indices, we obtain the following non- 
canonical Ricci tensor for this family of characteristic Killing tensors 



K 



■jk 



{ C*1212 + Cl313 + Cl414 Cl414 ~ Cl313 \ 

Ci414 — Ci313 C*1212 + Cl313 + C*1414 

2C 1313 + C 1212 

\ 2C UU + C 1212 ) 



After calculating the eigenvalues and corresponding eigenvectors of TZjk and applying 
the Gram-Schmidt orthonormalization procedure, we obtain an orthogonal matrix 



/ 


-V2 

2 








V2 
2 


\ 




V2 








V2 






2 


2 









-1 










V 








1 





/ 



A f 



which brings the Ricci tensor into canonical form 

Ttjk — diag(2C*i3i3 + C1212, 2C1313 + C1212, 2C1414 



+ C 



1212) 



2C4414 + C 



1212; 



26 



Therefore, we conclude that 



x = — ^j- cost cos u + ^- sin t sin v 
y = *y cos t cos u + ^ sin t sin v 
z = — cos t sin it 
w = sin i cos v. 

is a system of orthogonally separable coordinates for this Hamiltonian. 



6 Conclusion 

The results presented in this paper conclude an important project within the frame- 
work of a more general program of the development of the Hamilton- Jacobi theory 
of orthogonal separation of variables for natural Hamiltonians defined in spaces of 
constant curvature (see Table [3]). Having solved the equivalence problem, thus ex- 
tending the classical result by Olevsky, we have developed a general algorithm for 
solving the natural Hamiltonias defined in three-dimensional sphere via orthogonal 
separation of variables. In addition we give a simple and concise proof of the fact 
that the validity of the first and second TSN conditions imply the validity of the 
third. We have also derived a set of analogous algebraic conditions following from 
the vanishing Haanjes tensor which can be used to study and characterize algebraic 
and geometric properties of Killing two-tensors defined in spaces of constant, non-flat 
curvature. These conditions provide an alternative characterization of CKTs to the 
one derived in [21] and this paper based on the TSN criterion. 

The results presented here lay the groundwork for a project that concerns orthog- 
onal separation of variables afforded by characteristic Killing two-tensors defined in 
three-dimensional hyperbolic space which is the subject of a forthcoming paper. 
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Table 3: Solutions to the canonical forms and equivalence problems. 





Canonical forms problem 


Equivalence problem 


Euclidean space E 3 


Eisenhart, 1934 
Boyer et al, 1976 


Horwood et al, 2005 
Horwood, 2007 


Minkowski space M 3 


Horwood and 
McLenaghan, 2007 


Horwood et al, 2009 


Sphere 8 3 


Olevsky, 1950 


Cochran et al, 2010 



Appendix 

The following is a list of canonical forms for the six orthogonally separable coordinate 
systems of S 3 . Webs II, III, IV and VI are defined by K = aK\ + f3K 2 + 7C, where 
Ki, K 2 come from Eisenhart's equations, C is the Casimir tensor, and ct,/3,7 G R. 
Since it is possible to determine the presence of the Casimir tensor for webs I and 
V, and thus subtract it, these webs are defined by K = aKi + f3K 2 , where K x , K 2 
come from Eisenhart's equations and a,^6l. 

Rotational webs 

I. Spherical web 

K = Cl R 12 R l2 + 02(^13 R13 + #23 © R23) 
Rij = diag(ci + c 2 , c x + c 2 , 2c 2 , 0) 

II. Cylindrical web 

K = dR 12 R u + c 2 {R l3 R13 + Ru Ru + R23 R23 + #24 #24) + 

C3-R34 -R34 
Rij = diag(ci + 2c 2 , c x + 2c 2 , 2c 2 + c 3 , 2c 2 + c 3 ) 

III. Elliptic-cylindrical web of type 1 

K = Cl R 12 R 12 + c 2 (R 13 R 13 + R 23 #23) + c 3 (R u R 14 + R 24 #24)+ 

C4R34 -R 3 4 

Essential parameter: k 2 = — — 

c 4 - c 3 

Rij = diag(ci + c 2 + c 3 , c x + c 2 + c 3 , 2c 2 + c 4 , 2c 3 + c 4 ) 
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IV. Elliptic-cylindrical web of type 2 

K = CiR 12 #1 2 + C 2 (#i 3 #13 + #23 #23) + C 3 (#14 #14 + #24 #24) + 
C 4 # 34 #34 

Essential parameter: k 2 = — 

c 2 - c 3 

Rij = diag(ci + c 2 + c 3 , c x + c 2 + c 3 , 2c 2 + c 4 , 2c 3 + c 4 ) 

Translational web 

V. Sphere-elliptic web 

K = Ci# 12 #12 + C 2 #i3 #13 + C 3 # 23 # 23 

Essential parameter: k' 2 = — — 

ci - c 3 

Rij = diag(ci + c 2 , ci + c 3 , c 2 + c 3 , 0) 

Asymmetric web 

VI. Ellipsoidal web 

K = Ci#i 2 #i 2 + C 2 #13 #13 + C 3 #14 © #14 + C 4 #23 © #23 + C 5 #24 © #24 + 
C6#34 © #34 

where the q satisfy the constraint 

(c 3 + c 4 )(cic 6 - c 2 c 5 ) + (c 2 + c 5 )(c 3 c 4 - c 4 c 6 ) + (ci + c 6 )(c 2 c 5 - C 3 C 4 ) = 

Essential parameters: 
_ Ci(c2 - c 4 ) + c 6 (c 2 - c 3 ) - c 2 (c 3 + c 4 ) + 2c 3 c 4 

ci(c 2 - c 4 ) + c 4 (c 6 - c 2 ) + c 5 (c 4 - c 6 ) 
_ C 2 (C! - c 4 ) + ci(c 5 - c 4 ) - c 3 (ci + c 5 ) + 2c 3 c 4 
ci(c 2 - c 4 ) + c 4 (c 6 - c 2 ) + c 5 (c 4 - c 6 ) 
#jj = diag(ci + c 2 + c 3 , ci + c 4 + c 5 , c 2 + c 4 + c 6 , c 3 + c 5 + c 6 ) 
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